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A  novel  nonlinear  missile  guidance  law  against  maneuvering  targets  is  designed  based  on  the  principles 
of  partial  stability.  It  is  demonstrated  that  in  a  real  approach  which  is  adopted  with  actual  situations, 
each  state  of  the  guidance  system  must  have  a  special  behavior  and  asymptotic  stability  or  exponential 
stability  of  all  states  is  not  realistic.  Thus,  a  new  guidance  law  is  developed  based  on  the  partial  stability 
theorem  in  such  a  way  that  the  behaviors  of  states  in  the  closed-loop  system  are  in  conformity  with  a 
real  guidance  scenario  that  leads  to  collision.  The  performance  of  the  proposed  guidance  law  in  terms 
of  interception  time  and  control  effort  is  compared  with  the  sliding  mode  guidance  law  by  means  of 
numerical  simulations. 

©  2011  ISA.  Published  by  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Recently,  nonlinear  control  theories  have  been  used  in  the 
design  of  robust  guidance  laws.  Methods  such  as  Lyapunov-based 
nonlinear  guidance  laws  [1-3],  first-order  sliding  mode  guidance 
laws  [4-8]  and  nonlinear  guidance  laws  [9,10]  have  been 
obtained  from  Lyapunov  theorems  on  asymptotic  stability  or 
exponential  stability  of  all  states.  In  this  paper,  it  is  shown  that, 
in  practice,  such  a  behavior  is  not  realistic  for  all  the  states  of  a 
guidance  system. 

It  is  worth  noting  that  although  the  missile  guidance  problem 
can  be  treated  in  the  class  of  uncertain  nonlinear  systems,  it  has 
a  characteristic  property:  the  missile  should  intercept  a  highly 
maneuvering  target  within  a  finite  interception  time,  sometimes 
only  a  few  seconds.  In  this  paper,  a  guidance  law  will  be  designed 
which  guarantees  the  mentioned  characteristic.  The  present  work 
may  be  classified  in  the  category  of  robust  guidance  laws; 
however,  in  comparison  to  other  robust  guidance  laws,  it  follows 
a  completely  different  approach  towards  the  missile  guidance 
problem.  It  is  shown  that,  in  a  practical  approach  to  the  guidance 
problem,  each  state  has  a  specific  behavior  and  there  is  no  need 
for  asymptotic  convergence  of  all  the  states.  This  is  in  contrast 
to  conventional  methods  in  nonlinear  control  theory  which  try 
to  force  all  the  states  to  asymptotically  converge  to  the  origin 
(equilibrium  point).  Indeed,  this  paper  presents  a  new  approach  in 
nonlinear  guidance  law  design.  The  approach  is  advantageous  from 
a  practical  view-point  since  it  forces  the  states  to  behave  as  they 
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need  in  practice  and  leads  to  classification  of  the  state  variables 
of  the  guidance  system  dynamic  with  respect  to  their  required 
stability.  The  proposed  guidance  law  is  based  on  the  principle  of 
partial  stability ,  which  is  stability  with  respect  to  some  of  the  state 
variables. 

For  many  engineering  problems  like  inertial  navigation  sys¬ 
tems,  spacecraft  stabilization  via  gimbaled  gyroscopes  or  fly¬ 
wheels,  electromagnetic,  adaptive  stabilization,  etc.  [11-14],  the 
provision  of  partial  stability  is  necessary.  In  the  mentioned  appli¬ 
cations,  while  the  plant  may  be  unstable  in  the  standard  sense,  it 
is  partially  and  not  totally  asymptotically  stable.  The  effectiveness 
of  the  proposed  guidance  law  in  intercepting  a  highly  maneuver¬ 
ing  target  with  zero-miss-distance  and  finite  interception  time  is 
demonstrated  analytically  and  through  computer  simulation. 

2.  Partial  stability  analysis 

Partial  stability  is  defined  as  the  stability  of  a  dynamical  system 
with  respect  to  only  some  of  the  state  variables.  This  approach  is 
essential  for  many  engineering  fields  [14].  Consider  a  nonlinear 
dynamical  system: 

x=f(x),  x(t0)=x  o  (1) 

where  x  e  Rn  is  the  state  vector.  Let  vectors  Xi  and  x2  denote  the 
partitions  of  the  state  vector,  respectively.  Since  x  =  (x\,xl)T, 
Xi  e  Rni,  x2  e  R n2f  and  nA  +  n2  =  n,  the  nonlinear  system  (1) 
can  be  divided  into  two  parts: 

xi(t)  =  Fi(xi(t),  x2(t)),  xi(t0)  —  Xio 

x2(t)  =  F2(Xi(t),  x2(t)),  x2(t0)  =  X20  1  j 

where  x\  e  D  c  Rni ,  D  is  an  open  set  including  the  origin,  x2  e  Rn2 
andFi  :  DxR"2  ->  Rni  is  such  that  for  every  x2  e  Rn2,F^(0,x2)  =  0 
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and  Fi(.,  x2 )  is  locally  Lipschitz  in  x\.  Also,  F2  :  D  x  R"2  Fn 2  is 

such  that  for  every  Xi  e  D,  F2(xi,  .)  is  locally  Lipschitz  in  x2  and 
IXQ  =  [0,  tXq),  0  <  tXq  <  oo  is  the  maximal  solution  existence 
interval  (xi(t),  x2(t))  of  (2)  Vt  e  IXQ.  Under  these  conditions, 
the  existence  and  uniqueness  of  the  solution  is  ensured.  Now,  the 
stability  of  the  dynamical  system  (2)  with  respect  to  X\  can  be 
defined  as  follows  [14]. 

Definition  1.  (i)  The  nonlinear  system  (2)  is  Lyapunov  stable 
with  respect  to  X\  if  for  every  e  >  0  and  x20  e  Rn2,  there  exists 
5(£,  x2o)  >  0  such  that  ||xi0||  <  8  implies  ||xi(t)||  <  £  for  all 
t  >  0. 

(ii)  The  nonlinear  system  (2)  is  asymptotically  stable  with  respect 
to  Xi,  if  it  is  Lyapunov  stable  with  respect  to  X\  and  for  every 
x2o  G  Rn 2,  there  exists  8  =  8(x20)  >  0  such  that  ||xi0||  <  8 
implies  lim^ooXiCt)  =  0. 

In  order  to  analyze  partial  stability,  the  following  theorem  and  its 
corollary  are  taken  from  [14]. 


Theorem  1.  Nonlinear  dynamical  system  (2)  is  asymptotically  stable 
with  respect  to  x\  if  there  exist  a  continuously  differentiable  function 
V  :  D  x  Fn2  ->  R  and  class  I<  functions  a(.)  and  y(.),  such  that 


V(0,x2)  =  0,  x2eR"2 

a(lkill)  <  V(x,,x2),  (xi,x2)  e  D  x  R"2 

dV  (xi,  x2)  9V(x,,x2) 

-Fi(xi,x2)H - — - F2(x i,x2)  <  — y (||xt ||), 


3xi 

(xi,  x2)  g  D  x  Rf2. 


dx2 


(3) 

(4) 


(5) 


Proof.  See  [14].  □ 


Corollary  1.  Consider  the  nonlinear  dynamical  system  (2).  If  there 
exist  a  positive  definite,  continuously  differentiable  function  V  :  D  -> 
R  and  a  class  I<  function  y  (.),  such  that 

9V(xi) 

— r— —  Fi (xi,  x2)  <  -y (||X!  ||)  (x1?x2)  e  D  x  R  2  (6) 

3xi 

then  the  equilibrium  point  of  the  nonlinear  dynamical  system  (2)  is 
asymptotically  stable  with  respect  to  x\. 

Proof.  See  [14].  □ 


3.  Application  of  partial  stability  to  the  missile  guidance 
problem 


3.1.  Problem  formulation 

In  this  section,  a  two  dimensional  air  to  air  engagement  is 
considered.  The  vehicles  are  modeled  as  point  masses  and  the 
relative  kinematics  equation  of  motion  between  the  missile  and 
the  target  is  derived  in  the  planner  line  of  sight  (LOS)  coordinate 
system.  The  center  of  the  planner  LOS  coordinate  system  is  located 
on  the  missile.  The  first  axis  (XL)  is  along  the  LOS  vector  (i.e.  the 
vector  connecting  the  missile  to  the  target  at  each  moment 
and  its  direction  is  toward  the  target).  The  second  axis  (YJ  is 
perpendicular  to  the  first  one.  The  components  of  the  missile  and 
target  speed  vectors  {vM  and  vT)  in  the  LOS  coordinate  system  are 
shown  in  Fig.  1. 

The  representation  of  the  LOS  vector  in  the  LOS  coordinate 
system  is  as  follows: 


Fig.  1.  LOS  coordinate  system. 


where  r  is  the  relative  distance  between  the  missile  and  the  target. 
Also,  the  LOS  vector  makes  an  angle  of  0  with  respect  to  the  inertia 
reference  line.  The  LOS  coordinate  system  rotates  with  respect  to 
the  inertia  reference  line  and  con  -  defined  below  -  denotes  the 
rotation  rate  of  the  LOS-frame  with  respect  to  the  inertia-frame. 
The  only  nonzero  element  of  this  vector  is  perpendicular  to  the 
plane  of  the  LOS  coordinate  system: 


COIL 


0 

0 

<9 


(8) 


The  elements  of  the  relative  speed  vector  and  the  relative 
acceleration  vector  in  the  LOS  coordinate  are 


L 


Vr 


L(Vt  —  Vm)  = 
L0t  —  aM)  = 


Ki  —  Vm  i 

M 

\_VT2  ~  Vm2 

M 

\aj\  —  czmi 

an] 

| _aT2  —  aM  2 

_a,2J 

(9) 

(10) 


Now,  by  differentiating  with  respect  to  time  and  use  of  the 
Coriolis  Theorem  [15]  the  relative  equations  of  motion  in  the  LOS 
coordinate  can  be  obtained  as  follows: 


Lvr  =  l(D/?)  =  L(DLr)  +  lcoil  x  Lr 


r 

0 

r 

r 

0 

+ 

0 

X 

0 

= 

rb 

_0_ 

b 

_0_ 

0 

where  Lvr  is  the  relative  speed  vector  represented  in  the  LOS 
coordinate.  D/  is  the  differential  operator  with  respect  to  the  inertia 
coordinate  and  DL  is  the  differential  operator  with  respect  to  the 
LOS  coordinate. 

Note.  Since  the  inner  product  in  the  Coriolis  Theorem  is  used,  the 
LOS  vector  and  its  derivation  are  shown  with  three  elements.  Of 
course,  the  third  element  is  always  equal  to  zero. 

Therefore,  the  elements  of  the  relative  speed  vector  in  the  LOS 
coordinate  are 


M  _  ki  —  l’mi! 
rb  J  [ut2  —  VM2  \ 

where  r  and  rb  are  the  radial  and  tangential  components  of 
the  relative  speed  vector,  respectively.  Now,  by  differentiating 
the  relative  speed  vector,  the  acceleration  elements  in  the  LOS 
coordinate  can  be  obtained: 


Urll 

_^r2J 


(12) 


Lar  =  L(Djvr)  =  L(DLvr)  +L  coIL  x  % 


r 

0 

r 

r  -  rb2 

rb  +  rb 

+ 

0 

X 

rb 

= 

rb  +  2  rb 

0 

b 

0 

0 

r  = 


(13) 
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As  a  result,  the  elements  of  relative  speed  in  the  LOS  coordinate  are 


r  r  -  re2  1 

fln 

—  aM  i 

ar\ 

[i r'6  +  2 r9  ~\  ~ 

aT2 

—  0M2 

ar2 

(14) 


The  above  equations  are  defined  as  engagement  equations.  By 
considering  [Vr  Ve]  =  [r  rd],  Eq.  (14)  can  be  rewritten  in  the 
following  way: 


r ' 

Vr 

"0  0“ 

"0  0" 

Vr 

= 

Vg/r 

- 

1  0 

u  + 

1  0 

Ve. 

-VrVg/r 

_0  1_ 

o  i_ 

where  x  =  [r,  Vr,  Ve]T  is  the  state  vector,  w  =  (w^,w2)T  = 
(an,  Qt2)t  and  u  =  (ui,  u2)T  =  (aM i,  aM 2V  are  the  acceleration 
vectors  of  the  target  and  missile,  respectively.  Here,  the  accelera¬ 
tion  of  the  target  is  assumed  as  an  external  bounded  disturbance 
and  only  this  bound  is  required  in  the  design  of  the  guidance  law 
and  the  accurate  measurement  of  target  acceleration  during  the 
maneuver  is  not  necessary.  The  reason  for  choosing  rO  instead  of 
6  as  the  third  state  variable  is  that  it  avoids  the  appearance  of  the 
term  1  /r  in  the  control  and  disturbance  coefficient  matrix.  Thus, 
such  a  selection  would  make  these  matrices  constant. 


Remark  1.  The  initial  conditions  of  the  terminal  phase  are  usually 
in  such  a  way  that  r0  >  0  and  Vr0  <  0  (as  assumed  in  the  present 
work).  This  means  that  the  target  is  in  front  of  the  missile  and  the 
missile  is  approaching  it. 


3.2.  Desirable  behaviors  for  each  state 

For  the  interception,  it  is  sufficient  that  r(t)  becomes  zero  at  an 
instance  (r(t/)  =  0  where  tf  is  the  interception  time)  and  there 
is  no  need  for  r(t)  to  asymptotically  converge  to  zero.  In  other 
words,  asymptotic  convergence  is  not  an  ideal  behavior  for  r(t). 
It  should  be  noted  that  asymptotic  convergence  behavior  means 
that  the  missile  initially  approaches  the  target  very  fast;  however, 
near  the  target,  the  relative  distance  reduces  so  slowly  that  the 
missile  touches  the  target  in  an  infinite  time.  It  is  evident  that  such 
a  behavior  is  not  a  desirable  behavior  in  missile  guidance  which  is 
intended  to  hit  the  target  with  a  non-zero  speed  in  a  finite  time.  For 
this  purpose,  it  is  sufficient  for  the  relative  radial  speed  between 
the  missile  and  the  target  to  satisfy  the  following  proposition. 

Proposition  1.  In  order  to  intercept  the  target  within  a  finite  time 
interval,  it  is  sufficient  to  have 

3ti  e  [t0t/)  s.t.  Vr{t )  <  -?  <  0  Vt  e  [tit/].  (16) 


Proof.  r(t)  is  positive  and  continuous.  Consequently,  to  reach  a 
zero  value,  i.e.  zero-miss-distance,  within  a  finite  time,  it  should 
decrease  from  its  value  at  ti(<t/)  down  to  a  zero  value  at  tf. 
Inequality  (16)  implies  r(t)  —  r(ti)  <  — f(t  — H),  which  shows  that 
r  (t)  reduces  within  the  time  interval  [t\  t ] .  Since  it  is  desirable  to 
have  r(tf)  =  0,  thus 


t/  < 


r(ti)  +  fti 
K 


(17) 


Hence,  the  expectation  is  that  Vr  satisfies  Proposition  1.  In  other 
words,  the  convergence  and  stability  of  this  state  are  not  under 
consideration.  □ 


Remark  2.  A  common  approach  in  relevant  papers  is  to  regulate 
Vr  to  a  negative  constant  c.  Although  this  approach  guarantees  the 
interception  of  the  non-maneuvering  target  within  a  finite  time, 
it  is  not  efficient  for  intercepting  a  highly  maneuvering  target  in 
an  acceptable  interception  time.  In  such  a  case,  to  improve  the 


r 


performance,  c  should  vary  with  time;  however,  to  determine  it, 
the  maneuvering  model  of  the  target  should  be  known  for  the 
missile,  which  is  not  always  practically  possible. 

For  the  third  state  variable,  Ve,  asymptotic  convergence 
behavior  to  zero  is  desirable.  This  results  in  equality  of  vM2  and 
vT2f  that  is  equivalent  to  stating  that  the  LOS  vectors  remain 
parallel  with  each  other  (stopping  in  LOS  rotation).  The  following 
discussion  will  clarify  the  reason  for  the  acceptability  of  such  a 
behavior  [16]. 

Consider  the  first  equation  of  engagement  equations  (14): 

r  =  r02  +  (an  -  aM i).  (18) 

Assume  0  is  temporarily  fixed  and  the  resulting  radial 
acceleration  components  are  negligible.  Then  the  solution  of  this 
differential  equation  is  as  follows: 

r  =  c\eai  +  c2e~at  (19) 

where  C\  —  0.5 (r0  +  Vr0/a),  a  =  \6\  and  r0,  Vr0  are  the  initial 
conditions  of  engagement.  Usually,  these  initial  conditions  are  such 
that  r0  >  0  and  Vr0  <  0  (refer  to  Remark  1).  This  means  that  the 
target  is  in  front  of  the  missile  and  the  missile  is  approaching  it. 
Based  on  these  initial  values,  three  cases  might  occur  (in  the  third 
quarter  of  the  r  —  Vr  plane). 

-  If  Ci  <  0  or  r0  <  —  Vr0/a,  then  the  first  term  in  expression  (19) 
is  negative  which  results  in  reduction  of  the  relative  distance 
(r)  with  time.  Therefore,  the  missile  will  intercept  the  target. 
So,  this  area  is  called  the  “interception  region”. 

-  Ifci  >  Oorr0  >  —Vr0 /a,  then  the  relative  distance  increases 
with  time  and  the  missile  does  not  intercept  the  target.  This  area 
is  called  the  “non-interception  region”. 

-  If  Ci  =  Oorr0  =  — VrQ  / ot,  then  there  is  a  linear  relation  between 
r  and  Vr  at  all  moments.  This  line  is  the  boundary  between  the 
interception  and  non-interception  regions. 

These  regions  are  shown  in  Fig.  2.  As  observed,  it  is  obvious  that 
a  smaller  \6\  will  result  in  a  larger  interception  region.  Thus,  it  is 
desirable  that  6  converges  to  zero  (quickly)  and  stays  on  it.  As  a 
result,  asymptotic  stability  behavior  is  desirable  for  convergence  of 
6.  Also,  the  appropriate  behavior  for  Ve  is  the  same.  This  is  because 
stopping  in  LOS  rotation  leads  to  a  zero  value  for  Ve. 

Therefore,  the  state  vector  may  be  separated  into  X\  =  Ve  and 
x2  =  [r  Vr]  where  asymptotic  stability  behavior  only  for  X\  is 
desirable.  By  modeling  the  system  (14)  in  x\  —  x2  coordinates,  the 
following  can  be  obtained: 


*2 


-VrVe 

r 


—  u2  w2 


Vr 

V6 

- U\  +  ICi 

r 


(20) 


It  is  evident  that  only  the  second  component  of  the  input  vector 
appears  in  the  X\ -equation.  Therefore,  the  asymptotic  convergence 
of  Ve  may  be  achieved  by  u2.  In  the  meantime,  the  appropriate 
behavior  of  x2  could  be  obtained  by  u\. 
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3.3.  Guidance  law  design 

First,  consider  Eq.  (15)  where  w  =  0.  This  is  equivalent  to 
a  non-maneuvering  target.  Taking  the  Lyapunov  function  to  be 
V(x\)  =  0.5*1  —  0-5  V^2,  the  time  derivative  of  V(*i)  in  the  line 
of  the  system’s  trajectory  is 


Assume  |uq|  <  y) i,  and  take  iq  =  —r\\  sgn(Vr).  Since  Vr  is  sup¬ 
posed  to  be  negative,  thus  sgn(Vr)  =  — 1  and  iq  =  r\\.  By  inserting 
U\new  in  the  dynamic  equation  (15),  one  has: 

Vr  =  o  Vr  —  rj\  +  uq.  (31) 

Thus, 


VQci)  =  Vo 


( 


-VrVe 

r 


Take 


u2 


-VrVo 

r 


+  NVe. 


(21)  Vr(t)  =  Vr0eat  +  f  (— y/i  +  uq)ea(t  T)dr 
Jo 


(22) 


<  ^roear  - 

<  vr0eat. 


■  f  Thea«-r) dr  +  [‘me^dr 

Jo  Jo 


(32) 


Thus,  the  derivative  of  V(x\)  will  satisfy  the  condition  V(x\)  < 
—y  (||*i  ||)  where  y  (||*i  ||)  =  N  V2.  Therefore,  according  to  Corol¬ 
lary  1,  asymptotic  convergence  towards  zero  is  achieved  for  X\. 
Also,  by  choosing 


V2 

U\  = - oVr  G  >  0, 

r 

one  has  Vr  =  oVr,  and  consequently 


(23) 


Vr(t)  =  Vr0eat  (24) 

where  Vr0  <  0  (according  to  Remark  1).  As  a  result,  Vr(t )  <  Vr0 
<  0.  Therefore,  Proposition  1  is  satisfied. 


Remark  3.  Clearly,  higher  values  of  a  cause  a  shorter  interception 
time;  however,  its  adjustment  should  be  made  with  respect  to  the 
physical  limitations.  It  is  worth  noting  that  o  adjustment  can  also 
be  made  in  an  adaptive  manner  and  it  may  decrease  as  r  decreases. 
For  example,  it  can  decrease  linearly  as  follows: 


G(t)  = 


r0 


(To  + 


°f 


(25) 


where  a0  and  Of  denote  the  initial  and  final  values  of  a, 
respectively. 

Now,  assume  w  /  0.  In  this  case,  the  target  may  have  any  ar¬ 
bitrary  maneuver  with  a  bounded  acceleration.  Additional  control 
components,  v\  and  v2,  may  be  designed  in  such  a  way  as  to  make 
the  guidance  law  robust  with  respect  to  w.  By  taking 


Choosing  — f  =  Vr0  <  0  results  in  Vr(t)  <  — f  <  0  for 
t  e  [0  t/],  for  which  Proposition  1  is  satisfied. 


Remark  4.  Since  discontinuous  controllers  suffer  from  chattering, 
one  way  to  alleviate  this  problem  is  to  consider  an  approximation 
of  the  signum  function  by  a  saturation  function  with  a  high 
slope  (i). 

Consequently,  the  guidance  law  (33)  guarantees  interception  of 
the  maneuvering  target  within  a  finite  interception  and  zero-miss 
distance: 


Vr 

new(x)  =  —  -  oVr  +  m  or.  *7 1  >  0 

v0vr  fv6\ 

U2new(X)  ~ - - - ^  NVq  +  f]2Sat  l  —  \ 


(33) 

N,  s,  r\2  >  0. 


3.4.  Computer  simulations 

Numerical  simulations  are  performed  to  illustrate  the  effective¬ 
ness  of  the  designed  nonlinear  guidance  law. 

The  initial  conditions  of  engagement  are  specified  as 

r (0)  =  5000  m, 

Vr(0)  =  -300  m/s,  (34) 

Ve(0)  =  100/s 

and  the  designed  guidance  law  (33)  is  compared  with  the 
sliding  mode  guidance  law  presented  by  the  following  guidance 
command  [8]: 


u2new(x )  =  U2(x)  +  v2(x) 

VrVa 

=  — —+NV0  +  v2(x) 
r 


one  has 


(26) 


Ui  =  —  +  rj\  sat 
r 

u2  —  —  (N  +  1) - h  r]2  sat 

r 


(35) 


V(xi)  =  -NV 'l  -  Vg(v2  -  w2)  (27) 

where  the  last  term,  i.e.,  —Vq(v2  —  w2),  is  the  effect  of  the  control 
component  (u2)  and  disturbance  term  (w2).  Assume  \w2\  <  r\2\ 
therefore, 


-  V0(v2  -  w2)  <  -Vev2  +  |V0||u>2| 

<  -V0v2  +  rj2  \  Vq  I .  (28) 

By  choosing  v2  =  r\2  sgnCVgi),  one  has 

-  Vote  ~  w2)  <  -rf2 \Ve\  +  m \Ve\  =  0.  (29) 

Thus,  V(*i)  <  —y  (||*i  ||)  when  w2  is  present.  Now,  the  addi¬ 
tional  term,  iq,  could  be  designed  in  such  a  way  that  the  control 
law,  u^new(x)  =  Ui(*)  +  iq(*),  guarantees  the  specified  behavior 
for  *2  in  the  presence  of  uq. 

V2 

^neJx)  =  -*--trVr  +  VV  (30) 


The  sliding  mode  guidance  law  has  been  compared  with  some 
existing  guidance  laws  [8].  This  guidance  law  has  obtained  better 
performance  in  terms  of  interception  time  and  control  effort 
compared  to  other  existing  guidance  laws.  In  this  paper,  it  is  shown 
that  our  proposed  guidance  law  achieves  better  performance 
compared  to  the  sliding  mode  guidance  law.  In  both  cases,  a 
highly  maneuvering  target  with  the  following  acceleration  vector 
is  considered: 


aT(t)  = 


200  sin  —  t 
6 

200  sin  —t 
4  . 


(36) 


The  initial  conditions  in  the  inertia-frame  of  the  missile  and  the 
target  for  nonlinear  simulation  are  given  by  the  following. 


-  For  the  missile:  *M(0)  =  yM(0)  =  0  m. 

-  For  the  target:  *T(0)  =  5000  cos  00  m,yT(0)  =  5000  sin  00  m, 
vTx(0)  =  —200  m/s,  vTy( 0)  =  0  m/s. 
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Fig.  3.  (a)  Trajectories  of  the  missile  and  target,  (b)  Relative  radial  speed  (Vr). 


Fig.  4.  (a)  Radial  guidance  command  (u\ ).  (b)  Tangential  guidance  command  ( 112 ). 

Oq  is  set  to  be  7r/6.  The  initial  values  of  the  missile  speed 
components  can  be  obtained  based  on  the  initial  relative  speed 
components  (34)  and  the  initial  values  of  target  speeds. 

The  constant  N  is  chosen  as  4,  while  r\\  and  rj2  are  set  to  be  200. 
Also,  s  =  1,  =  0.1,  a  =  0.05  and  c  =  -800  are  selected. 


Fig.  3(a)  shows  the  trajectories  of  the  missile  and  the  target  and 
Fig.  3(b)  depicts  the  relative  radial  speed  (Vr).  Fig.  4(a)  and  (b) 
reflect  the  guidance  commands. 

Fig.  3(a)  illustrates  the  collision  points,  Cl  and  C2,  for  the 
proposed  guidance  law  and  the  sliding  mode  guidance  law, 
respectively.  The  interception  time  is  7.26  s  for  the  proposed 
guidance  law  and  8.38  s  for  the  sliding  mode  guidance  law.  Fig.  4(a) 
demonstrates  that  the  amount  of  change  of  radial  control  effort  for 
the  proposed  law  is  less  than  for  the  sliding  mode  guidance  law. 
Also  Fig.  4(b)  shows  that  the  proposed  guidance  law  requires  less 
tangential  control  effort  than  the  sliding  mode  guidance  law. 

4.  Conclusion 

In  this  paper,  a  new  approach  to  the  missile  guidance  problem 
was  presented.  It  became  clear  that  in  a  successful  missile 
guidance  scenario,  which  leads  to  target  interception,  the  desirable 
behaviors  of  the  state  variables  are  different  from  each  other 
and  asymptotic  convergence  behavior  is  not  ideal  for  all  the 
state  variables.  Therefore,  based  on  the  partial  stability  theorem, 
a  new  robust  missile  guidance  law  was  developed  and  its 
effectiveness  in  finite  time  interception  of  maneuvering  targets 
was  analytically  shown.  Also,  the  proposed  guidance  law  provided 
better  performance  in  comparison  with  the  sliding  mode  guidance 
law. 
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